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Abstract. Considering the Casimir effect due to phononic excitations of a 
weakly interacting dilute BEC, we derive a re- normalized expression for the zero 
temperature Casimir energy Eq of a BEC confined to a parallel plate geometry 
with periodic boundary conditions. Our expression is formally equivalent to the 
free energy of a bosonic field at finite temperature, with a nontrivial density of 
modes that we compute analytically. As a function of the interaction strength, £c 
smoothly describes the transition from the weakly interacting Bogoliubov regime 
to the non-interacting ideal BEC. For the weakly interacting case, Sc reduces 
to leading order to the Casimir energy due to zero-point fluctuations of massless 
phonon modes. In the limit of an ideal Bose gas, our result correctly describes 
the Casimir energy going to zero. 



PACS numbers: 67.10.Fj, 03.75.Hh, 42.50.Lc 
1. Introduction 

The Casimir effect is a consequence of the distorted vacuum fluctuation spectrum of 
quantized fields in bounded domains or spaces with non-trivial topologies [1,2]. In 
Casimir's original calculation, the system under consideration is the electromagnetic 
vacuum [3]. Imposing Dirichlet boundary conditions along one spatial direction by 
confining the system between two (idealized) parallel plates causes a change in the 
(infinite) vacuum energy-density. The variation of the vacuum energy-density with 
respect to the plate separation is called Casimir pressure, and, after renormalization, 
yields a finite expression for an attractive interaction energy per unit area between 
the plates. The electromagnetic Casimir force, caused by quantum fluctuations of the 
electromagnetic vacuum, is varying as he/ with the plate separation L and the speed 
of light c. It has been measured in a number of experiments using various experimental 
settings [4, sec. 6]. The comparison between quantum vacuum experiments like 
these and the predictions of different theoretical models provides possibility to test 
fundamental physics (like higher dimensions or additional interactions), in much the 
same way as accelerator experiments in high-energy physics do at the other end of the 
energy scale [5,6]. 

We consider here, instead of the electromagnetic vacuum, a weakly interacting 
Bose-Einstein-Condensate (BEC) at zero temperature, and expect, in a similar 
manner, the quantum fluctuations on top of the ground state of the BEC to give 
rise to observable Casimir forces: Within the Bogoliubov approximation, the excited 
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states of a BEC can be treated as quasi-particles characterized by the dispersion 
relation 



that behaves Hncar for small momenta, with the 'soimd velocity' c = h/{2mC,) being 
inverse to the healing length C,. In (1.1), the wave-number 1/^ characterizes the 
transition between the linear (phonon) and the quadratic (free-particle) regimes [7] . Q 
is also related to the s-wave scattering length a of the atoms and to the BEC density 
n via C, = l/(4y^7ma). For small momenta, the quasi-particles (phonons) propagate in 
the same way as the massless electromagnetic field, except for the propagation velocity 
being different. Hence, the zero temperature quantum fluctuations in a spatially 
confined BEC can be expected to result in an observable Casimir force. 

Different scenarios for Casimir forces in BECs have been analyzed by previous 
work: for the parallel plate geometry, an asymptotic expansion of the Casimir force has 
been calculated in [8] , the small expansion parameter being the ratio between healing 
length and plate separation. In the leading order, it reproduces exactly the same 
he/ L'^ behavior as in the electromagnetic vacuum. The next order corrections scale 
with the ratio C/L- Replacing the perfectly reflecting plates by impurities embedded 
in a quantum liquid, Casimir forces between these impurities have been calculated 
in [9, 10] as a function of the impurity-liquid coupling. If the impurities are realized by 
atoms which, in a certain internal state, interact with the atoms of the quantum 
liquid through s-wave scattering, the Casimir interaction should be detectable as 
a shift of spectral lines that depends on the distance between the impurities [10]. 
In the limit of an infinitely strong impurity-liquid coupling, the result for the (one- 
dimensional) parallel plate scenario was recovered [9]; for a weak coupling, however, the 
interaction between the impurities vanishes exponentially with the impurity separation 
on a scale set by the healing length. For the ideal (i.e. non-interacting) Bose gas, 
it was found that there is no Casimir force at all between impurities of arbitrary 
interaction strength, including the idealized parallel plate scenario [9, appendix C]. 
This is consistent with the quite general method of [11], which is mapping (polynomial) 
dispersion relations to Casimir forces in the parallel plate geometry: this method shows 
that media with quadratic dispersion relations, and hence the ideal Bose gas, do not 
give rise to any zero-temperature Casimir forces. A non-vanishing Casimir force in the 
ideal Bose gas can arise due to thermal fluctuations, as calculated in [12,13]. All these 
forces are small but finite observable quantities, which if cxpc^rimentally confirmed — 
would provide direct evidence of the quantum fluctuations in weakly interacting BECs. 

The system under consideration in the present article is a homogeneous, weakly- 
interacting dilute BEC at zero temperature, confined to a parallel plate geometry with 
periodic boundary conditions in one of the three spatial dimensions. The condition 
of diluteness can be formulated as n\a\'^ <^ 1, where a is the s-wave scattering length 
and n = N/V is the density, N being the total particle number and V the volume 
of the gas. The quantity n\a\^ is usually called the gas parameter. The perturbative 
calculation of typical properties of such a BEC, like the ground state energy or the 
depletion of the condensate, is then essentially an expansion in the gas parameter [7] . 

For this system, we will give a re-normalized expression for the Casimir energy- 
density per unit area. Our expression has the form of an integral over a 'density of 
modes' p{x) times the Bose distribution function: 



E{k) = hc^/k'^{l + k^^C) 



(1.1) 




(1.2) 
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where p{x) has a simple analytic form (soc (3.11)), and correctly describes the 
vanishing of the Casimir force in the limit of the interaction strength going to zero. 
The possibility to express the zero temperature Casimir energy in the above form, 
resembling the density of states of a bosonic system at finite temperature, is connected 
to a topological analogy between our parallel plate scenario and finite temperature 
field theory: In the parallel plate geometry, one spatial coordinate of the field is 
subject to periodic boundary conditions, while in finite temperature field theory, the 
imaginary time coordinate is subject to a similar periodicity condition. This analogy 
has been pointed out some time ago in [14]; it does not carry over, however, to Casimir 
calculations for non-linear dispersion relations. The periodic boundary conditions 
for the BEC have mainly been chosen because they make the relation to the finite 
temperature case particularly evident. In experiments, periodic boundary conditions 
can be realized in toroidal traps, but they also appear in optical lattices. If the 
boundaries arc taken as real physical plates, the perfect mirror scenario (i.e. Dirichlct 
boundary conditions as discussed in [ 1 1 , 1 2 , 1 5] ) is closer to an experimentally realizable 
situation. When Dirichlet boundary conditions are imposed on the fluctuations on top 
of the BEC ground state, the expansion of £c for moderate interaction shows in the 
leading term again the same behaviour as a massless scalar field propagating at the 
speed of sound c, but with a different numerical prefactor. 

The article is organized as follows: In section 2, we briefly recall some well- 
known perturbative expressions for the free energy, chemical potential and ground 
state energy of a dilute, weakly interacting BEC at zero temperature. Working from 
the formula for the free energy, we derive in section 3 the expression for the Casimir 
energy. In the limit (/L <^ 1, the function p reproduces the power-series expression for 
the Casimir energy derived in [8], as is shown in section 3.1. In section 3.2, we relate 
the above mentioned analogy to the finite temperature case in some detail, as this topic 
does not seem to get overly much attention in the recent literature. Finally, section 3.3 
deals with the behavior of the function p in the limit of the effective coupling constant 
between the particles going to zero. Here, our expression £c correctly describes the 
vanishing of the Casimir force. We comment on the failure of the large-distance 
expansion to describe the non-interacting Bose gas. The case of Dirichlet boundary 
conditions is briefly discussed in the Appendix. 

2. Free energy in a weakly interacting, dilute BEC at T = 

In the one-loop approximation, the (unrenormalized) free energy density of a weakly 
interacting BEC at T = (in three spatial dimensions) is given by 



Here and in the following, we work in units where 2m = H = ks = 'i-- In the above 

expression, /i is the chemical potential (which, in section 3.1, will be connected to 
the speed of sound in the medium), and g is the effective coupling constant that can, 
for low energy collisions in a dilute medium, be identified with the s-wave scattering 



with the mean-field and one- loop contributions [16]: 



(2.1) 




(2.2) 



(2.3) 
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length o: g = Sira. The term A^J^ represents the one- loop counter-term needed to 

render J^i finite. If the integral in (2.3) is regularized with a momentum cut-ofF in 
the ultraviolet, the linear, cubic and quintic UV- divergences in J^i can be absorbed by 
re-normalizing g, ji and the vacuum energy respectively. If dimensional regularization 
is used in (2.3). no counter-term is needed at this level of perturbation theory. 

After renormalization, the chemical potential ji can be obtained from (2.1) by 
inverting 

n(.) = -f, (2.4) 

which yields 

IJL{n) = gn{l + 0{Vna^)) . (2.5) 

Reinserting /x into (2.1) will then reproduce the classical result for the leading quantum 
corrections to the ground state energy density £ of a hard-sphere Bose gas, which was 
first derived by Lee, Huang and Yang in [17]: 



|l + (na^)'/' + 0{naHn[na^])^ (2.6) 



3. The Casimir contribution as an integral over a density-function 

Now, in order to describe a BEC between a pair of parallel plates separated by a finite 
distance L (with periodic boundary conditions), we have to quantize the momentum 
component perpendicular to the plates: 

27r 

k'^ k'^+ujl, ujn^—n, neZ (3.1) 

Correspondingly, the momentum integration perpendicular to the plates is replaced 
by a discrete sum: 




The area of the plates is taken to be L1L2 with Li,I/2 3> L, so the system now 
inhabits the volume V = L1L2L. The mean field contribution J^o (see (2.2)) depends 
only trivially on the new boundary conditions, with the volume V entering through n 
in (2.5). 

The Casimir energy £c of the BEC is related to the free energy at the one-loop 
level. We are interested in its change per unit area that is due to the introduction of 
the boundary conditions: 

= LJ'il^^y + Sc (3.3) 

where the first term gives the (one-loop contribution to) the free energy in a 
homogeneous system. After applying the substitutions (3.1) and (3.2) to (2.3), we 
are left with the following expression for J^i, now describing the leading quantum 
corrections to the free energy of a BEC confined between parallel plates: 

^i = -=2 E ^n(^k)'+n^][M{kr+n^], (3.4) 

where we have used the abbreviations 

M{k,L,C)=L^/WTT/e , (3-5) 
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(a) Area of integration in (3.7). (b) 'Density of states' function p{x). 

Figure 1. 

Fig. 1(a): The upper curve shows the function M(fc, L, Q from (3.5), the lower one 
is Lk = Lfe/(27r). The two shaded regions make up the total area of integration 
in (3.7). 

Fig. 1(b): The exact 'density of states' function p(x) from (3.11) (solid line) and 
its small-argument expansion (3.12)) (broken line). 



with L = L/{2tt) and ( — l/y^. The summation over n can be converted into two 
integrals by using the Abel-Plana formula in the form (see [18]) 

f: fin) = r dx fix) + I /(O) + ^ Tdx ^ ^'"^'J '^''"'^ (3.6) 
n=0 -^0 ^ "^0 e -1 

Application of (3.6) splits (3.4) into the two terms written in (3.3), as the /(O)-tcrm 
cancels out. The first term, L!Fi\y_y, is divergent, but of the same form as the 
one-loop free energy (2.3) itself. Hence, it can be re-normalized as described above. 

The term £c in (3.3), an energy per area, describes the effects induced by 
restricting the periodicity in one spatial dimension to be much smaller than the 
remaining ones, i.e., confining the system to a volume V with one 'short' side L. 
£c is convergent for any finite L, strictly negative, and goes to zero for L — > oo: 

__2 [ d-k [x--{Lkf]"' [M{kf-x-]"' 

The integration domain is sketched in figure 1(a) (the shaded regions). This domain as 
well as the form of the integrand in (3.7) are due to the branch-points of the integrand 
in (3.4): In order to stay clear of the branch-cuts, we have evaluated the last term in 
the Abel-Plana formula (3.6) along an integration contour slightly to the right of the 
imaginary axis. Inserting the function / given by (3.4), 



f{x) = \I{x^+l' k^)ix^+AP) , (3.8) 
the numerator in the last term of (3.6) evaluates to 

i{f{ix + 0) ~ f(-ix + 0)) = -2^(.T2-L^fc2)(M2-a;2), Lk<x<M, (3.9) 

being zero everywhere else along the contour of integration (see [1, sec. 2.2]). By 
changing the order of integration in (3.7), the fc-integral can be performed, and we 
find the expression (1.2) for the Casimir energy: 
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Figure 2. Casimir energy per unit area as a function of the plate separation 
L = L/(2tt), normalized to the healing length Recall that C, = ^^1/ (gnrn) 
with g the effective interaction constant and n the BEC density. The Casimir 
energy has been normalized to its value £sc{L) = —{TT^/90)flc/L^ for a masslcss 
scalar field propagating with a velocity c = h/{2mC,); this limit is approached 

at large distance ( ). The broken curve ( ) shows the asymptotic 

expansion for L/^ ^ 1 (,3.1,3), while the dotted curve ( — • — ) shows the opposite 
limit L/C, — > (non-interacting Bose gas), of eqn. (3.23). The full black curve is 
obtained by numerically evaluating the integral in (3.10) with the mode density 
p{x) of (3.11). It smoothly describes the dependence of Sq on the interaction 
strength in the regime where both of the asymptotic expansions diverge. 



In the form of (3.10), £c is expressed as an integral over a 'density of states' (DOS) 
for a bosonic system, as setting the integration variable x = /3a;/(27r) will reproduce 
the Boltzmann factor in the denominator. 

The 'mode density' p{x) is obtained by integrating separately over the lower 
triangular region in figure 1(a) (dark gray) and the 'hyperbolic tail' (light gray). In 
terms of the dimcnsionless variable r] = (C/L) x, we have 



Above X = L/C, the function p{x) changes into a constant independent of x (see the 
solid line in figure 1(b)). The Casimir energy £c as a function of the normalized 
distance L/C is shown in figure 2, obtained by numerically integrating (3.10). 

To provide a cross-check for the above results, we will in the following section show 
that the function p{x) reproduces the asymptotic expansion for the Casimir energy 
given in [8, 19], which is valid for small values of the parameter C/L- In section 3.3 we 
then show that the formula for p(x) in (3.11) also yields correct results in the opposite 
limit of the ideal BEC, i.e. C ~^ oo. 

3.1. Expansion at large distance and moderate interactions 

In this section, we will assume the plate separation to be much greater than the healing 
length of the BEC, i.e., the ratio C/L can be treated as a small parameter. Note 




for < a; < L/C 



for X > L/C 



(3.11) 
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that this hmit cannot describe a strongly interacting Bose gas in the proper sense, 
since this would lead to the break-down of the one-loop approximation at the basis of 
our approach. 

As can be seen in figure 1(b), the kink in the function p{x) happens at the large 
value X = L/C, 3> 1 in our Umit. In evaluating the Casimir energy with (3.10), large 
values oix, and hence the behavior of p{x) after the kink, get exponentially suppressed 
by the denominator. Thus we can approximate the function p{x) by expanding its 
small- a; part as a power-series around a; = 0: 

The integrals in the expression (3.10) for the Casimir energy can be performed 
explicitly.^ Upon setting c := l/Q (the speed of sound in the medium), we recover the 
result derived in [8]: 

The leading term £sc = — 7r^c/(90-L^) is, as noted in [8,15], equal to the Casimir 
energy (per unit area) of a massless scalar field with propagation velocity c confined 
between two parallel plates with periodic boundary conditions [20]. Its presence is a 
manifestation of the Goldstone-theorem, the long wavelength part of the Bogoliubov 
spectrum representing the gap-less Goldstone modes [21]. The next-to-leading term in 
(3.13) is in [8] referred to as the Bogoliubov correction due to the non-linearity of the 
dispersion. Indeed, if we express the Bogoliubov dispersion relation (1.1) as a power 
series around A; = and put the first few terms into the expression for the Casimir 
force derived in [11, eqn. (16)], the leading term in the same manner reproduces the 
result for the scalar field (in one dimension), followed by terms that are smaller in 
magnitude and of opposite sign. As can be seen in figure 2, the contribution of these 
corrections is getting smaller as the ratio between plate separation and healing length 
increases, leaving only the dominant scalar term £sc depicted by the horizontal dashed 
line in figure 2. 

When finally calculating the Casimir pressure from (3.13), one has to consider 
that, as the number of particles is held constant, the derivative of the speed of sound 
with respect to L is not zero, as discussed in [8]. 



3. 2. A formal analogy to finite temperature systems 

As already mentioned, our expression (3.10) for £c at T = formally resembles the 
DOS for a bosonic system at finite temperature. This can be understood by recalling 
that a finite temperature system can be described in imaginary time, combined with 
periodic boundary conditions with period (3. The same topology is realized in the 
parallel plate system at T = 0, when one spatial dimension is subject to periodic 
boundary conditions with period L ( [14], see also [22]). 

X We use the identities 

^oo ^2r-i ^ r(2r)C(2r) 
Jo e'^^t - 1 ~ (27r)2'- 
in terms of the Gamma- and Zeta-functions. 
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Let us briefly re-phrase the argument of [14] within our notation. The canonical 
partition function Z for a system at temperature T = 1/(3 in D spatial dimensions is 
usually expressed as the path integral 

r0 



j>d[^ exp dT j d^xC{^) 



(3.14) 



where C is the Lagrangian (a scalar functional of the field $), and the field is 

constrained in such a way that $(x, 0) = $(x, /3). If wc take C = (c^/2) (?^<I> 9''$, 
the path- integral in (3.14) can, after Fourier expansion of the field, be evaluated to 
yield [23,24] 

\n[Z] = -V-Y,jj^H^l + ^m^ ^n = jn, u;{k)=c\k\. (3.15) 

Here, LUn with n = 0, ±1, ±2, . . . arc the Matsubara-frcqucncics due to the periodicity 
condition in (3.14), and u>k with the continuous parameter k is the dispersion relation of 
the massless scalar field. Note the similarity to (3.1), where we had periodic boundary 
conditions not in imaginary time but in one spatial dimension. 

The sum over n in (3.15) is usually evaluated by multiplying with a factor 
|/3 cot(|/3a;), which has poles of residue 1 at w = 27rn//3, and then integrating over 
a contour in the complex cj-plane which includes all the poles [24, sec. 3.4]. This 
technique is actually the same that is used in [18] to prove the Abel-Plana formula 
(3.6). The well known result is 



ln[Z] = \n[Z] 



/3->oo 

where we have already subtracted the zero-point fluctuations. After integrating by 
parts and employing the integral-identity from the footnote on page 7, we obtain (fi^j 
is the volume of the unit sphere in D dimensions) 



(2,)' 

Note that, upon setting k — (27r//3c) x, the above expression has the same form as £c 
in (3.10), with p{x) being proportional to x^. Now, let the system inhabit a volume 
V = L1L2L. With D = 3 we get for the free energy per unit area 

= r3i8i 

LiL2 90 /3(/3c)3 ■ ^ ' ' 

Tentatively exchanging /3c with L in (3.18) will reproduce the Casimir energy for a 
massless scalar fleld, i.e. the first term in (3.13) [14]: 

L1L2 90 ^ ^ 

So, we have seen that the zero temperature Casimir energy of a massless scalar 
field confined between two parallel plates can be obtained by a simple exchange of 
variables, once we know the thermal contribution to the free energy density of that 
field. Unfortunately, this simple mapping does not carry over to fields characterized by 
nonlinear dispersion relations: The temperature dependent part of the free energy for 
a BEC is (in the Bogoliubov approximation) still given by (3.16) with the dispersion 
now being co{k) = ck{k'^ + [16]. But a simple interchange of f3c with L 

in (3.16) will no longer yield £c, because the discretized Matsubara frequencies a;„ 
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(that become a discretized momentum) always enter in (3.15) in the same way as a 
spatial momentum component enters in a linear dispersion relation. It thus seems 
very difScult to mimic the fully nonlinear behavior of the dispersion relation. 



3.3. The non-interacting limit 

The transition from the weakly interacting BEC (here described at the one-loop level 
only) to the ideal BEC should be accomplished by letting the effective coupling g = Swa 
(see section 2) go to zero, corresponding to C ^ cxd. The Casimir energy is expected 
to vanish in this limit, as shown in [9, 11]. 

The series expansion (3.12) was constructed for (/L <C 1, which is a physically 
reasonable assumption for finite L and a weak but finite interaction. But with the 
effective coupling strength g ^ 0, the chemical potential /u in (2.5) will vanish, too, 
and the healing length ( will diverge. As L is kept finite, the non-interacting BEC is 
hence described by the limit L/( 0, which is the opposite to the case considered 
in section 3.1. The asymptotic form of (3.12) for the mode density p{x) — as well 
as the Casimir energy in (3.13) — diverges in the limit of zero interaction strength. 
Our calculation of the exact mode density p{x) suggests that this divergence is due 
to a branch point in the complex x-plane that moves towards a; = and makes the 
power series expansion behind (3.13) break down. We show here that the exact mode 
density (3.11) leads to a Casimir energy that smoothly vanishes with the interaction 
strength (see figure 2). 

To examine the asymptotic behavior of £c for L/( <C 1, we again start from 
(3.10) and (3.11), separately treating the behavior of p{x) to the left and the right of 
the kink at L/<^: 

_ F'^ dxpjx) 1 dx 

= £i+£2 (3.21) 

As the upper limit of the integration in £i is going to zero, we can replace p{x) by the 
first term in the expansion of (3.12), and Taylor-expand the denominator about x = 0. 
After integrating over this expansion, £i will yield terms of 0({L/Q'^y Integrating 
£2, from the lower border of the integral we get a contribution 

Altogether, we find for the behavior of £c in the non-interacting limit 

= - T l^e^^^'^/^^ - 1] } + 0{(Lm ■ (3.23) 

Note that the leading order for C ^ 00 at fixed L goes like C~^liiC. Conversely, at 
fixed C) a logarithmic divergence remains for L 0. The logarithmic term changes 
sign for L/(^ > In[2] w 0.7 and, for large values of L/C, the above expression diverges, 
just as the expansion (3.13) does for small values of L/C (see figure 2). Hence, (3.23) 
and (3.13) provide two asymptotic expansions to £c for opposite limits, while the 
exact formula is given by (3.10) integrated over (3.11). 
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4. Summary 

Starting from the free energy in a weakly-interacting dilute BEC, we derived a re- 
normalized expression for the 'phononic' Casimir energy of the BEC confined at zero 
temperature to a parallel plate geometry with periodic boimdary conditions. Our 
formula for the Casimir energy (per unit plate area), (3.10), has the form of an integral 
over a density function p times the Bose distribution. The function p is given by a 
rather simple analytic expression in (3.11). In section 3.1, wc provided a cross-check 
for our result by showing that a series expansion of p in the parameter L/C, ^ 1 
reproduces the asymptotic series for the Casimir energy derived in [8]. There, the 
Euler-MacLaurin formula was used to extract the long wavelength behavior out of 
the UV-divergent sum over all Bogoliubov modes satisfying the boundary conditions. 
This approach fails to reproduce the non- interacting limit. 

As pointed out in [8,15], the Casimir energy of the weakly interacting BEC is, due 
to the linear dispersion of its low lying excitations, in the leading order determined 
by a term analogous to the Casimir energy of a massless scalar field propagating with 
vacuum velocity c = 1/^. Our result displays this behavior in the regime of the weakly 
interacting BEC where the plate separation is much larger than the healing length, as 
can be seen in figure 2. In addition, for C — > oo with L kept finite (the non- interacting 
limit), our result (3.11) correctly describes the Casimir energy going to zero and 
displays the Casimir energy as a smoothly varying function of the interaction strength 
in the intermediate range. The subtleties of the asymptotic expansions illustrate 
the rich physical content behind the nonlinear dispersion relation of the Bogoliubov 
vacuum. 
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Appendix: Dirichlet boundary-conditions 

Here we discuss how our results for the Casimir energy get modified with Dirichlet 
boundary conditions instead of periodic ones: taking the wavenumbers a;„ in (3.1) as 

ujn = nn/L with n running from 1 to infinity, application of the Abel-Plana formula 
generates an additional divergent surface-term Esurf in (3.3); 



(The same term with opposite sign will appear for von Neumann boundary conditions, 
similar to [13] where the ideal gas at finite temperature is considered.) As Esurf docs 
not depend on L, it does not affect the Casimir energy. Our main formula for Ec-, 
(3.10), is modified by an overall factor i, while in the mode density p{x, L) (see (3.11)) 
the distance argument picks up an additional factor of two: 




(A.1) 




(A.2) 
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The asymptotic expansion for large plate separation and moderate interactions, as 
derived in section 3.1, now yields (compare with (3.13)): 

Again, the leading term is equal to the known result for the Casimir energy of 
a massless scalar field with propagation velocity c = l/( confined between two 
parallel ideal mirrors (see [15, eqn. (38)]). For the non-interacting limit, as treated in 
section 3.3, we get the asymptotic expression (compare with (3.23)) 

This goes to zero as for periodic boundary conditions, with slightly different numerical 
factors. 

The divergent surface energy (A.l) does not fit into the renormalization scheme 
for the bulk part of the free energy J^i as described after (2.3). For a full treatment of 
the free energy of a BEC between perfect mirrors, it is not sufficient to change only the 
wavenumbers u)n in the Bogoliubov excitations, because both the ground-state wave 
function and the excitations have to be zero on the boundaries. A detailed analysis 
will be reported elsewhere. 
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